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We have measured the beta-decay branching ratio for the transition from 21 Na to the first excited 
state of 21 Ne. A recently published test of the standard model, which was based on a measurement 
of the f3-v correlation in the decay of 21 Na, depended on this branching ratio. However, until now 
only relatively imprecise (and, in some cases, contradictory) values existed for it. Our new result, 
4.74(4)%, reduces but does not remove the reported discrepancy with the standard model. 
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PACS numbers: 27.30.+t, 23.40.-s 



I. INTRODUCTION 



A recent publication by Scielzo et al. pj reported a mea- 
surement of the j3-v angular correlation coefficient, ap v , 
for the /3-decay transition between 21 Na and the ground 
state of its mirror, 21 Ne. The authors compare their re- 
sult with the standard-model prediction for apv, with a 
view to testing for scalar or tensor currents, the presence 
of which would signal the need for an extension of the 
standard model. Although they found a significant dis- 
crepancy - the measured value, ap v = 0.524(9), disagrees 
with the standard-model prediction of 0.558 - they stop 
short of claiming a fundamental disagreement with the 
standard model. 

Scielzo et al. pj offer two alternative explanations that 
would have to be eliminated before their result could 
begin to raise questions about the need for an exten- 
sion to the standard model. One is that some 21 Na2 
dimers formed by cold photoassociation could also have 
been present in their trap, thus distorting the result; they 
themselves propose to do further measurements to test 
that possibility. The other is that the branching-ratio 
value they used for j3 decay to the first excited state of 
21 Ne might not be correct. Because Scielzo's measure- 
ment could not distinguish between positrons from the 
two predominant /3-decay branches from 21 Na (see Fig. 

the adopted branching ratio for the (3 transition to the 
first excited state not only affects their data analysis but 
also helps determine the theoretical prediction for ap v it- 
self, since the axial- vector component of the ground-state 
branch can only be determined from its ft value, which 
also depends on the branching ratio. This branching ra- 
tio is a key component of their standard-model test, yet 
the five published values [2j, |3j, Lft la, |6j are between 25 and 
45 years old, are quite inconsistent with one another and 
range from 2.2(3) to 5.1(2)%. To remedy this problem, 
we report here a new measurement of the ground-state 
branching ratio, for which we quote ±0.8% relative pre- 
cision, five times better than the best precision claimed 



in any previous measurement. 



II. EXPERIMENT 

We produced 22.5-s 21 Na using a 28^-MeV 22 Ne beam 
from the Texas A&M K500 superconducting cyclotron to 
initiate the 1 H( 22 Ne, 2n) 21 Na reaction on a LN 2 -cooled 
hydrogen gas target. The ejectiles from the reaction were 
fully stripped and, after passing through the MARS spec- 



jpec 

trometer 0, produced a 21 Na secondary beam of >99% 
purity at the extraction slits in the MARS focal plane. 
This beam, containing ^3 x 10 5 atoms/s at 24.4A MeV, 
then exited the vacuum system through a 50-/iin-thick 
Kapton window, passed successively through a 0.3-mm- 
thick BC-404 scintillator and a stack of aluminum de- 
graders, finally stopping in the 76-//m-thick aluminized 
mylar tape of a tape transport system. Since the few im- 
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FIG. 1: /3-decay scheme for 21 Na 
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FIG. 2: Spectrum of /3-delayed 7 rays observed in coincidence 
with positrons following the decay of 21 Na. It includes about 
half of the total data collected. The peaks attributable to 
21 Na are marked with their energy in keV; the sum peak is 
identified by its components. The small unmarked peak at 
682 keV is caused by summing of one 511-keV 7 ray with the 
back-scattered 7 ray (171 keV) from the second 511-keV 7 
ray. 



purities remaining in the beam had ranges different from 
that of 21 Na, most were not collected on the tape; resid- 
ual collected impurities were concluded to be less than 
0.1% of the 21 Na content. 

In a typical measurement, we collected 21 Na on the 
tape for a few seconds, then interrupted the beam and 
triggered the tape-transport system to move the sample 
in 180 ms to a shielded counting station located 90 cm 
away, where the sample was positioned between a 1-mm- 
thick BC404 scintillator to detect /3 + particles, and a 70% 
HPGe detector for 7 rays. Two timing modes were used: 
in one, the collection and detection periods were 3 and 
30 s, respectively; in the other, they were 6 and 60 s. In 
both cases, after the detection period was complete, the 
cycle was repeated and, in all, some 3,200 cycles were 
completed over a span of 32 hours. 

Time-tagged /3-7 coincidence data were stored event by 
event. The (3 and 7-ray energies, the coincidence time be- 
tween them, and the time of the event after the beginning 
of the cycle were all recorded, as was the total number of 
/3-singles events for each cycle. The same discriminator 
signal used for scaling was also used in establishing the 
/3-7 coincidences. 

Essential to our experimental method is the precise 
absolute efficiency of the 7-ray detector, which was po- 
sitioned 15 cm from the collected sample. We have 
meticulously calibrated our HPGe detector at this dis- 
tance over a five-year period using, in total, 13 individual 
sources from 10 different radionuclides: 48 Cr, 60 Co, 88 Y, 

108m Ag) 109 Cd; 120m Sb; 133^ 134^ 137 Cg and ISOmjjf 

Two of the 60 Co sources were specially prepared by the 
Physikalisch-Technische Bundesanstalt [j| with activities 
certified to 0.06%. The details of our calibration pro- 



cedures, which include both source measurements and 
Monte Carlo calculations, have been published elsewhere 
H G3 E3 The absolute efficiency of our detector is 
known to 0.2% in the energy range from 50 to 1400 keV, 
and to 0.4% from 1400 keV to 3.5 MeV. 

The absolute efficiency of the j3 + detector, which was 
located 1.5 cm from the collected sample, is not required 
for our measurement but its dependence on energy is of 
some importance (see section llTT)l . We have explored the 
efficiency of this detector via measurements and Monte 
Carlo calculations, and its dependence on /3 + energy is 
now reasonably well understood [T^ | . 

A typical 7-ray spectrum recorded in coincidence with 
betas is presented in Fig. [21 Apart from the annihilation 
radiation, the only significant peak in the spectrum is 
the 351-keV 7 ray from the first excited state in 21 Ne. In 
3,200 total cycles we recorded more than 8 x 10 counts 
in this peak. 

It was important to our later analysis that we estab- 
lish the contribution of room background both to the 
/3-7 coincidence spectrum and to the /3-detector singles 
rate. For this purpose, we recorded data with the cy- 
clotron beam on but with a thick degrader inserted just 
upstream from the tape; everything was thus identical 
to a normal measurement except that no 21 Na was im- 
planted in the tape. Both the coincidence and singles 
rates were observed to drop to 0.04% of the rate observed 
when 21 Na was correctly implanted. Room background 
was thus effectively negligible in our analysis. 



III. RESULTS 

The /3-decay scheme of 21 Na is shown in Fig. ^ The 
only branch in addition to those populating the ground 
and first excited states is known to be very weak [6j and 
can be ignored in our analysis. In that case, the branch- 
ing ratio, Ri, for population of the first excited state can 
be determined from the measured intensity ratio of the 
351-keV 7 ray relative to the total number of 21 Na de- 
cays. Thus, we obtain R\ from the following relationship: 



Ri 



N 10 



k, 



(1) 



where N~p is the number of 351-keV 7 rays observed in 
coincidence with betas; N/3 is the number of (singles) be- 
tas observed; e 7 is the efficiency of the HPGe detector for 
351-keV 7 rays; and k is a factor (~ 1) that accounts for 
small experimental corrections that will be enumerated 
in what follows. Note that the efficiency of the beta de- 
tector does not appear in Eq. ^ although its dependence 
on [3 + energy will be seen to play a minor role in the 
evaluation of k. 

Before determining the ratio NJ/3/N/3 from our data, 
we eliminated those cycles in which the collected source 
was not positioned exactly between the /3 and 7 detec- 
tors. Although the tape-transport system is quite con- 
sistent in placing the collected source within ±3 mm of 
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the designated counting location, it is a mechanical sys- 
tem, and occasionally larger deviations occur. For each 
cycle we recorded not only the total number of positrons 
detected but also the total number of 21 Na ions that 
emerged from the MARS spectrometer, as detected by 
the scintillator located immediately in front of the alu- 
minum degraders. The ratio of the former to the latter 
is a very sensitive measure of how well the source is po- 
sitioned with respect to the /3 detector. In analyzing 
the data, we rejected the results from any cycle with an 
anomalous (low) ratio. Under these conditions, we ob- 
tained the result N y p/Np = 2.378(13) x 1(T 4 . 

As stated in section [HI the absolute efficiency, e 7 , of 
our detector at 15 cm is known to ±0.2%. However, 
this applies to a highly controlled situation in which the 
source-to-detector distance can be measured by microm- 
eter to a small fraction of a millimeter. With the fast 
tape-transport delivery system, we cannot be assured of 
reproducibility at the same level of precision. Taking 
±0.5 mm to be our actual uncertainty in position un- 
der experimental conditions, we add an uncertainty of 
±0.6% to the detector efficiency in quadrature with the 
basic ±0.2% uncertainty. For the 351-keV 7 ray, this 
leads to e 7 = 5.12(3) x 10~ 3 , the value we insert in Eq. 

m 

Although the ratio N 1 p/Np and e 7 are the predom- 
inant experimental quantities required to evaluate the 
branching ratio, it is the correction factor k that holds 
the key to our achieving high precision. In fact, k is re- 
ally a product of four separate corrections, ki...k±. We 
will deal with each individually. 

Random coincidences (fci). — Since the time between 
each coincident (3 and 7 ray was recorded event by event, 
we could project out the time spectrum corresponding to 
the 351-keV 7 ray. In that spectrum, the prompt coinci- 
dence peak stood prominently above the flat random dis- 
tribution, allowing us clearly to distinguish the relative 
contributions of real and random coincidences. The cor- 
rection factor required to account for the random contri- 
bution to the /3-coincident 351-keV 7-ray peak was thus 
determined to be k\ — 0.9884(10). Naturally, this cor- 
rection accounts not only for random coincidences among 
21 Na (3 and 7 rays but also for random coincidences be- 
tween 21 Na betas and any 7 rays originating from room 
background. 

Real- coincidence summing (fe). — Since each 351-keV 
7 ray from the decay of the first excited state in 21 Ne 
is accompanied by a positron from the 21 Na /3 + -decay 
branch that populated the state, there is a significant 
probability that a 351-keV 7 ray and 511-keV annihila- 
tion radiation will reach our HPGe detector simultane- 
ously and be recorded as a single 7 ray with the combined 
energy of both. Any summing of this kind will rob events 
from the 351-keV photopeak. Our first step in account- 
ing for the resultant loss was to obtain the area of the 
observed 862-keV (511±351) sum peak. Since losses from 
the 351-keV photopeak result not just from its summing 
with the 511-keV photopeak but also with the latter's 



Compton scattered radiation, as a second step we mul- 
tiplied the sum-peak area by the known "total-to-peak" 
ratio for our detector at 511 keV (see Fig. 11 in reference 
10]). Finally, this result for losses was increased by 4% to 
account for annihilation in flight, which leads to 351-keV 
peak summing with annihilation radiation of different en- 
ergies, and by another 2.5% to account for summing with 
positrons backscattered from the plastic scintillator. The 
total loss due to real-coincidence summing was thus de- 
termined to be 1.78%: i.e. k 2 = 1.0178(17). 

Dead time {k^). — Equation ^ depends upon N^p and 
Np being recorded for identical times. In our experiment 
they were, of course, gated on and off together, but dur- 
ing the counting period the circuit dead time for N 7 p, 
which was limited by the relatively slow electronics used 
for 7-ray counting, was much greater than that for Np, 
which was simply scaled. We determined the dead time 
associated with N 1 p from the total rate in the HPGe de- 
tector during the counting period and from the known 
processing time (32 /xs) for each coincident event. The 
scaler dead time per event was only 100 ns but the total 
rate in the scaler was much higher than the HPGe rate; 
nevertheless the dead time associated with Np turned out 
to be smaller by a factor of three than that associated 
with coincidence events. The overall correction factor is 
k 3 = 1.0018(1). 

Beta-detector response function (fei). — The correction 
factor associated with the /3-detector response function is 
given by 

ki = tPtotai _ 0.956ft, + 0.05e ft 
eft % 

where ep and ep ± are the detector efficiencies for the f3 
transitions to the ground and first excited states respec- 
tively. If the detector response function were completely 
independent of energy, then this correction factor would 
be unity. In fact, though, the efficiency does change 
slightly with energy. We have studied this effect using 
measurements with sources - 90 Sr, 133 Ba, 137 Cs and 207 Bi 
- aided by Monte Carlo calculations [ijj]. Including the 
effects of our low-energy electronic threshold, we deter- 
mine that ep 1 /ep a = 0.987, which leads to the correction 
factor fc 4 = 1.0129(13). 



TABLE I: Error budget for the measured branching ratio Ri 



Origin of uncertainty 


% uncertainty 


Experimental ratio, N 1 p /N/3 


0.52 


HPGe detector efficiency 


0.20 


Source-detector distance 


0.60 


Random coincidences 


0.11 


Real-coincidence summing 


0.17 


Dead time 


0.01 


/3-detector efficiency vs energy 


0.13 


Total uncertainty on Ri 


0.85 
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Multiplying ki through k^ we determine the correction 
factor in Eq. [J to be k = 1.0208(24). When combined 
in Eq. ^ with the other factors already discussed, this 
yields the final result for the branching ratio to the first 
excited state in 21 Na: 

Rx = 0.0474(4). (3) 

The complete error budget corresponding to our quoted 
±0.85% uncertainty is given in Tabled 



IV. ANALYSIS 

Our measured branching-ratio value is compared with 
previous measurements in Table [n] All previous experi- 
ments determined the branching ratio from a comparison 
of the area of the 351-keV peak to that of the annihila- 
tion radiation. This method has the advantage that only 
relative detector efficiencies are required, but it has three 
serious disadvantages: i) contaminant activities may well 
make an unknown contribution to the annihilation radi- 
ation; ii) most positrons do not annihilate at the source 
position, where the 7 rays originate, so the relative detec- 
tion efficiencies cannot be simply determined from cali- 
bration sources; and iii) the significant effect (~5%) of 
positron annihilation in flight is a first-order correction 
that must be calculated and corrected for. All previ- 
ous measurements except possibly reference Q were sus- 
ceptible topotential contaminants; only the last three 
references El 013 mention accounting for a spatially dis- 
tributed source of 511-keV radiation; and only the last 
two [5J, |6j appear to have taken account of annihilation 
in flight. 

Given the age of the previous measurements and 
the potential hazards associated with their experimen- 
tal method - not to mention their mutual inconsistency 
- we choose not to average our result with them but in- 
stead to use our present result alone in extracting the 
properties of the 21 Na /3-decay scheme. 

Since there are only two significant /3-decay branches 
from 21 Na - to the ground and first excited states of the 
daughter - with R\ determined, the branching ratio to 
the ground state, Rq, follows directly from it: 





#0 = 0.9526(4), 


(4) 


TABLE II: Measurements of the branching 


ratio _Ri 


Date 


Reference 


Result(%) 


1960 
1963 
1974 
1977 
1980 


Talbert & Stewart [2] 
Arnell & Wernbom [3] 
Alburger [4] 

Azuelos, Kitching & Ramavataram [5J 
Wilson, Kavanagh & Mann [6] 


2.2(3) 
2.3(2) 
5.1(2) 
4.2(2) 
4.97(16) 


2006 


This measurement 


4.74(4) 



where this result is actually determined to a precision of 
0.04%. We now proceed from this value for R to obtain 
the ft value for this transition, the relative contributions 
of axial-vector and vector components, and ultimately 
the standard-model expectation for its (3-v angular cor- 
relation coefficient. 

In deriving the ft value for the ground-state mir- 
ror transition, we take the half-life of 21 Na to be t-yi^ 
= 22.49(4) s and its total decay energy to be Qec = 
3547.6(7) keV. The former is the average of two mutu- 
ally consistent results 0, and the latter is the value 
quoted in the 2003 Atomic Mass Evaluation 0] where 
it was obtained from a single 20 Ne(p,7) 21 Na measure- 
ment made in 1969 and then revised by Audi et al. 
|13| to take account of more up-to-date calibration en- 
ergies. With the calculated electron-capture probability 
for the ground-state transition being 0.00095, the average 
half life, when combined with our branching ratio value 
from Eq. 0] yields a partial half-life for the transition of 
23.63(4) s. 

Next we compute the value of / from the Qec value 
following methods similar to those we used in the anal- 
ysis of super allowed decay; these are described in the 
Appendix to reference 15]. To make an "exact" cal- 
culation that includes, for example, the effects of weak 
magnetism and other induced corrections we need a shell- 
model calculation of the appropriate nuclear matrix ele- 
ments. For this we used an (s, <i)-shell model space and 
the universal (s, <f)-shell effective interaction of Wilden- 
thal [l6|. This interaction has been demonstrated 0] 
to reproduce energy spectra and Gamow- Teller matrix 
elements in this mass region providing that the axial- 
vector coupling constant is quenched. In our calculation, 
we fine-tuned the amount of quenching to reproduce our 
experimental data 1 . 

For a mirror transition like this one, which includes 
both vector and axial-vector components, the / value 
calculated for the vector part of the weak interaction, / v , 
is slightly different from the value calculated for the axial- 
vector part, f A . In the allowed approximation it is always 
assumed that f v = f A = f but, where high precision is 
sought, a more exact calculation is required. The results 
we obtain, / v = 170.974 and f A = 174.157, are nearly 
2% different from one another, principally as a result of 
the influence of weak magnetism on the shape correction 
factor of the axial- vector component. In quoting the ft 
value for the mirror ground-state transition, we make the 
(arbitrary) choice to use / v , with the result that 

f v t = 4040(9)s. (5) 

Like any other ft value, this result can be related to 
vector and axial- vector coupling constants, and to the 



We adjusted the quenching so that it reproduced our measured 
value of A (see Eq. IIUI . This corresponded to Ga — Gv and 
is essentially the same result that Brown and Wildenthal |17| 
established for the shell as a whole. 
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matrix elements pertaining to the specific transition. To 
do so with the precision required for a standard-model 
test requires that radiative and charge-dependent cor- 
rections be incorporated. The expression we use is the 
following: 



f v Gl(l) 2 (l + 6 NS -5 c )(l + AZ) 
+/ A G 2 » 2 (1 - 6 A )(l + A£) 1 (1 + S' R )t = K (6) 



where K/{hcf = (8120.271 ± 0.012) x 10~ 10 GeV~ 4 s; 
G v and G A are the vector and axial-vector coupling con- 
stants for nuclear weak decay; and (1) and (a) are the 
Fermi (vector) and Gamow- Teller (axial-vector) matrix 
elements, respectively, for the ground-state transition. 
For this particular transition between T—\ states, (1) = 1. 
The transition-dependent radiative correction terms, 5' R 
and 5ns, and the isospin-symmetry-brea king correction, 
5c, all have their conventional definitions |15| but, in the 
present context of a mixed vector and axial-vector tran- 
sition, we note that 5' R is the same for both components 
while 5c and 5ns on ly pertain to the vector component. 
The latter two terms have their equivalents that must be 
applied to the axial-vector component but we subsume 
them into a term we call 5a'- as it turns out, we will not 
have to calculate a value for 5a- Finally, the transition- 
independent radiative correction also takes on different 
values for the vector and axial- vector components, A^ 
and A*; but neither will have to be calculated. 
Rearranging Eq. |BJ we obtain the result: 



f v t(l + 5' R )(l + 5 NS -5c) 



K 



where 



A = 



G, 



W(l-^) 1/2 (1 



G 2 (1 + A X ) [1 + A 2 £] 

A A)l/2 



(7) 



G v (l + <5 J vs-'5c) 1/2 (l + Av)i/2- 



Here A is the ratio of axial- vector to vector components in 
the transition. A further simplification in this equation 
can be achieved by our implementing the results from 
superallowed + — » + beta decays, which provide an 
experimental determination of the product G 2 (l + A^). 
The average corrected Ft value from these decays [l5| is 
related to the vector coupling constant via the relation- 
ship: 



Tt 



K 



2G 2 (1 + Av)- 



(8) 



Since it is the term A that we need to extract from 
experiment in order to calculate the (3-v angular correla- 
tion coefficient, we now re-express Eq. \7\ in the following 
form: 



A 



k 

f, 



ITt 



f v t(l + S' R )(1 + 5 NS - 6c) 



1 



(9) 



We have calculated the three remaining correction terms 
using the same methods as were described in reference 



[3, the results being S' R = 1.492(15)%, 6 C = 0.268(16)% 
and 5 NS = -0.065(20)%. We then adopt the value, Tt 
— 3072.7(8), which is the average result extracted from 
superallowed + — * + beta decays when the correction 
terms are calculated by the same methods as those used 
here (see Eq. 11 in reference 15] ). Thus we finally obtain 



A = 0.7033(24) 



(10) 



for the ground-state mirror transition. 

Based on this result for A we have computed the beta- 
neutrino correlation coefficient exactly, following the for- 
malism of Behrens-Biihring |19| . These authors write the 
electron-neutrino correlation u>(9, W) as: 



u(6, W) = ^D{k,W)P k {cos( 



(11) 



where W is the electron energy (in rest-mass units), 9 
is the angle between the emitted electron and neutrino 
directions, and Pk are Legendre polynomials. The sum is 
over k — 0, 1, 2. The coefficients D(k, W) are expressed 
fully by Behrens and Buhring |l9j, from which it can be 
seen that D(0,W) is exactly equal to 1.0, and D(2,W) 
is small. The term D(l, W) relates to the beta-neutrino 
angular correlation coefficient, ap v , via the expression 



D(1,W) = a 0vP /W, 



(12) 



where p = \/W 2 
we compute 



1. For the exact expression for ap u 



a^ = (D(l,W)W/p), 



(13) 



where (..) signifies an average over the beta spectrum. It 
should be noted that this exact evaluation of ap^ yields 
a result that is about 1% different from the approximate 
expression that is often used: viz. 



approx 



= (l-A 2 /3)/(l + A 2 ). 



(14) 



The exact expression in Eq. 1131 differs from this approx- 
imate one by the inclusion of energy dependence as well 
as weak magnetism and other small effects. Our final 
computed result for the exact f3-v angular correlation co- 
efficient based on our new experimental result for A is 



a Pu = 0.553(2) 



(15) 



This can now stand as the "standard-model prediction" 
for apv, against which the measured angular-correlation 
coefficient can be compared. Our new value is 0.9% lower 
than the one originally used by Scielzo et al. 0. 



V. CONCLUSIONS 

As noted in the Introduction, the value of the branch- 
ing ratio affects not only the standard-model prediction 
for ap v (see Eq. I15|) but also the analysis by Scielzo tt 
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al. [l| of their measurement of that coefficient. With the 
excited-state branching ratio taken to be 5.02(13)%, they 
applied a correction of +6.81(18)% to their result. Since 
this correction scales with the branching ratio pfj|. our 
new value for the latter leads to a new correction factor 
of +6.44(5)%. This downward shift of 0.4% is actually 
rather small compared to the overall uncertainties quoted 
by Scielzo et al., and their value for ap v as obtained from 
21 Ne 1+ only changes from 0.524(9) to 0.523(9). 

As a result of our new measurement we have improved 
- and lowered slightly - the standard-model prediction 
of the (3-v angular correlation coefficient for the mirror 
transition from 21 Na. This new prediction still leaves 
the Scielzo et al. experimental result Q] in disagreement 
with the prediction. However, the authors themselves 
expressed concern about the possible presence of 21 Na2 
dimers in their trapped samples; this would have caused 
a dependence of their result on the trapped-atom pop- 
ulation and could easily reconcile their result with the 
standard model. With a precise branching ratio now de- 



termined, an investigation of the actual make-up of the 
trapped-atom samples in the Scielzo et al. experiment is 
essential if the 21 Na result is to become a real test of the 
standard model. 

Finally, we draw attention to the fact that the 
"standard-model prediction" for ap v depends on the half- 
life of 21 Na and its /3-decay Q value through the f v t value 
for the ground-state transition (see Eq. [S] and the pro- 
ceeding paragraphs). The half- life has only been mea- 
sured twice 0, |5| - in experiments that did not obtain 
branching ratios in agreement with our current result - 
and the Q value comes from a single 35-year-old (p,7) 
measurement 0] originally based on long-outdated cal- 
ibration energies. Clearly, both these results could be 
improved significantly by modern measurements. 

The authors would like to thank N.D. Scielzo for 
helpful discussions. This work was supported by the 
U.S. Department of Energy under Grant No. DE-FG03- 
93ER40773 and by the Robert A. Welch Foundation un- 
der Grant No. A-1397. 
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